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基于 修正 后 矩阵 分 解 的 最 优 协 方差 DOA 估计 
外 ” 敬 ， 马 黎 文 ， 申 “ 东 ， 蒋 占 军 ， 李 琴 然 


(兰州 交通 大 学 电子 与 信息 工程 学 院 , 兰州 730070) 


摘 要 : 针对 传统 来 波 方 向 〈direction-of-arrival，DOA) 估计 在 信号 相干 、 低 信 品 比 与 噪声 非 均 匀 环 境 下 性 能 差 的 
问题 ， 基 于 修正 后 的 算 阵 分 解 ， 提 出 一 种 利用 凸 优化 的 协 方差 矩阵 最 优 DOA 估计 方法 。 修 正 后 的 矩 阵 分 解 方法 ， 
解 相 干 的 同时 克服 了 孔径 损失 ; 然后 ， 利 用 凸 优化 ， 重 构 出 无 噪声 的 协 方差 矩阵 ; 最 后 ， 利 用 最 小 化 搜索 计算 出 
DOA。 仿真 结果 表明 ,所 提 算 法 与 矩阵 分 解 (matrix decomposition, MD) 算法 、 基 于 工 范 数 的 奇异 值 分 解 (Li-norm 
singular vector decomposition，Li-SVD) 算法 以 及 基于 空间 平滑 的 协 方差 秩 最 小 化 估计 (spatial smoothing based 
covariance rank minimization，SS-CRM) 算法 比较 ， 能 更 好 地 抑制 非 均匀 噪声 ， 且 在 低 信 噪 比 条 件 下 ， 依 然 性 能 让 
好 。 
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Optimal covariance matrix based on modified matrix decomposition algorithm 


Di Jing, Ma Liwen, Shen Dong, Jiang Zhanjun, Li Cuiran 
(School of Electronics & Information Engineering, Lanzhou Jiaotong University, Lanzhou 730070, China) 


Abstract: To solve the problems that the traditional direction-of-arrival(DOA) estimations have poor performance when 
processing coherent signals in the cases of low signal-to-noise ratio and non-uniform noise . This paper proposed a DOA 
method based on the modified matrix decomposition method, with the best covariance matrix by the convex optimization. 
The modified matrix decomposition method can deal with the extraction of the coherent sources, while, overcomes the 
aperture loss. Furthermore, the method reconstrtucted the noise free covariance matrix by the convex optimization. Finally, 
using minimization search to calculate DOA. The simulation results show that, comparing to the matrix decomposition (MD) 
algorithm, Li-norm singular vector decomposition(Li-SVD) algorithm and spatial Smoothing based covariance rank 
minimization (SS-CRM) algorithm, the method suppress the non-uniform noise well and have good performance in low 
signal noise ratio (SNR). 
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0 ”外 言 空间 并 不 准确 0223 ,其 信号 估计 的 性 能 并 不 能 得 到 保障 。 为 

一 了 进一步 提高 DOA 算法 的 性 能 ，Malioutov 等 人 提出 了 基于 
在 无 线 通 信 中 ， 由 于 多 径 信 号 相干 的 影响 ， 信 和 号 协 方差 Z1 范 数 的 奇异 值 分 解 (Li-norm singular vector decomposition， 
CS 和 矩阵 会 出 现 秩 缺 现象 ， 使 信号 子 空间 的 维 数 小 于 信号 源 数 ， Li-SVD ) 稀 玻 重 构 的 DOA 估计 0 , 较 之 传统 的 子 空间 方法 ， 
将 信号 子 空间 “扩散 ”到 噪声 子 空间 由 ， 导 致 一 些 相 干 信号 稀疏 重 构 方法 能 够 同时 处 理 非 相 干 及 相干 信号 ， 具 有 更 高 的 
源 的 导向 矢量 与 噪声 子 空 间 不 完全 正 交 ， 给 正确 估计 信和 号 源 角度 分 辨 力 ， 在 少量 快 拍 数 以 及 低 信 品 比 条 件 下 有 更 高 的 估 
方向 造成 困难 中。 同时 ， 在 实际 应 用 中 ， 由 于 一 些 信道 的 硬 计 精 度 。 但 是 它 在 低 信 品 比 情况 下 ， 得 到 的 解 矢量 稀疏 性 较 
牛 差 异 与 非 期 望 干扰 ， 使 得 实际 的 噪声 模型 并 不 满足 高 斯 差 ， 空 间 谱 中 存在 较 多 伪 峰 ， 对 DOA 的 正确 


噪声 ， 直 接 导 致 传统 高 分 辨 率 的 算法 性 能 恶化 中 。 为 了 解决 a(0.)= 
这 个 问题 ， 针 对 相干 信号 的 来 波 方向 (DOA) 估计 算法 相继 Tr 27 ，， 2r ， .oT (1) 
a es de 训 Le -ji 生 sina ) i exp Jj 人 Easing || 
被 提出 ， 这 些 算法 通常 都 利用 阵列 导向 矢量 的 特殊 性 质 ， 对 L 和 和 
协 方 差 窍 阵 的 秩 进行 恢复 ， 从 而 达到 解 相 干 的 目的 上 5 。 则 阵列 输出 信号 可 写成 如 下 的 矩阵 : 

经 典 的 解 相 干 算法 有 空间 平滑 类 算法 [6-51, 和 矩阵 重 构 类 算 X=AS+G 
法 &10 等 。 和 矩阵 分 解 算 法 (matrix decomposition, MD) 属于 古 计 带 来 了 一 定 的 干扰 1 中 。 而 在 实际 应 用 中 ， 由 于 接受 
和 矩阵 重 构 类 算法 ， 虽 然 可 以 解决 信号 相关 性 的 问题 ， 但 是 ， 设备 的 性 能 差异 与 外 界 影响 ， 各 个 传感器 的 噪声 通常 都 是 不 
MD 算法 的 解 相 干 能 力 是 通过 减少 矩阵 的 维 数 来 换取 的 00。 相干 的 非 均 匀 噪 声 ， 会 对 传统 的 DOA 估计 算法 的 性 能 带 来 
在 低 信 噪 比 以 及 非 均匀 噪声 的 条 件 下 ， 信 和 号 协 方差 矩阵 的 刺 很 大 的 影响 9。 针对 该 问题 ， 文 献 [17] 提 出 一 种 基于 空间 平 


角 线 元 素 也 是 非 均匀 的 ， 所 以 由 特征 矢量 直接 构造 的 信号 滑 的 协 方差 秩 最 小 化 估计 (spatial smoothing based covariance 
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rank minimization，SS-CRM) 算 法 , 虽然 可 以 有 效 的 降低 
的 影响 ， 然 而 ， 在 相 邻 信 
是 无 法 有 效 检测 出 信号 DOA。 
域 下 加 权 Li 稀 玻 导 


匀 品 声 对 信 
小 信 噪 比 过 低 的 情况 下 , 还 
文献 59 提出 的 二 阶 统计 量 
( weighted 


号 子 空间 


Sparse 


reconstruction 


所 教 ， 等 : 


DOA 


基于 修正 后 和 琵 阵 分 解 的 


FE 均 
号 源 角 度 差 较 


而 


EE 构 DOA 估计 
estimation 


algorithm, Li1 -MC-WLOSRSS ) 算法 虽然 对 非 均 匀 噪 声 有 很 好 


的 抑 杀 


能 力 ， 但 是 加 


MY 


也 带 
针对 上 述 问 题 


来 了 限制 。 


种 基于 


大 减少 奇异 值 分 解 计算 
过 凸 优 化 原理 ， 


， 本 文 将 会 
修正 后 抵 阵 分 解 的 最 优 

(optimal covariance matrix based 
decomposition, OC-SMD) 。 
平滑 技术 中 的 前 向 与 后 向 


权 算 法 也 带 来 了 


以 矩阵 分 解 算法 为 基础 ， 提 
办 方差 矩阵 DOA 估 


量 的 同时 提 


1 信 


OO 
1 2 


元 位 


标 信 号 源 ， 入 射 


da 了 且 


6 。 系 统 是 由 W 个 阵 元 构成 的 等 
分 别 为 0.4,…， 


间距 直线 阵 ， 各 
1 所 示 。 


M>K， 如 图 


阵 元 M-1 


阵 元 M 


图 1. 信 号 模型 


Fig.1 Signal model 


较 大 的 计算 量 ， 对 于 其 


[a 
LU 


对 


on the smoothed matrix 
修正 后 的 矩阵 分 解 就 是 利用 空 
平滑 矩阵 来 重 构 协 方差 矩阵 ， 在 大 
高 算法 的 稳定 性 。 然 后 通 
求 出 非 均匀 的 噪声 功率 ， 继 而 得 到 无 噪声 的 
协 方差 矩阵 。 最 后 进行 最 小 化 搜索 09 实 现 信号 的 DOA 估计 。 


号 模型 
假设 有 xX 个 远 场 窄带 


间 


角度 分 别 为 


阵 


Q= diag[q9,,92,..., 


最 优 协 方差 DOA 估计 


元 上 的 非 均匀 噪声 功率 。 


式 (5) 可 知 ， 可 以 通过 
bb 信和 号 的 DOA。 但 是 当 信 
亏 ， 所 以 特征 


值 分 解 从 而 求 吕 
接收 信号 的 协 方差 矩阵 利 
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9,] 是 非 均匀 噪声 协 方差 矩阵 ，9, 指 第 m 个 阵 


对 接收 信 


号 协 方差 矩阵 进行 特征 


值 少 于 信号 源 数 


特征 


， 所 以 无 法 实 


写 源 为 相干 
值 分 解 之 后 得 到 的 的 
岗 DOA 估计 。 


言 号 时 ， 


2 OA 
i 


2.1 
假设 


修正 后 的 协 方 差 矩 阵 分 解 算法 
一 个 均匀 线 阵 有 M 个 阵 元 ， 分 成 P(P=M-4+D 个 


子 阵 ,每 个 子 阵 具 有 9 个 阵 元 。 另 设 有 天 个 信号 源 入 射 ， 信 和 号 
协 方差 矩阵 为 式 (5)， 则 第 个子 阵 的 接收 信号 协 方差 矩阵 为 
R=R(n:n+g-l)。 则 前 向 平滑 P 次 的 数据 协 方差 算 阵 为 : 
] 2 
R= (0) 
则 后 向 平滑 P 次 的 数据 协 方差 矩阵 可 表示 为 
R=JRY (7) 
J 为 9x4 的 对 角 和 矩阵 。 
则 将 数据 协 方差 矩阵 Rw 作 如 下 改进 : 
Rs =[R.,R,] (8) 
对 上 述 协 方 差 矩阵 进行 奇异 值 分 解 (SVD) 就 可 得 到 修正 


后 的 信号 子 空间 与 
本 文 提 到 的 修 
decomposition, SMD) 算 法 。 
但 是 ,划分 子 阵 的 做 法 
SMD 在 奇异 值 分 


二 
VT 污 - 屋 


占 计 分 辨 紊 的 降低 。 
办 方差 的 特征 分 解 所 带 来 的 信号 子 空间 泄漏 2" 


虽然 能 够 降低 信 


DOA 估计 性 外 


EE 的 大 幅 下 降 。 


多 匀 品 


2.2 基于 SMD 的 最 优 协 方差 矩阵 DOA 估计 


为 J 


做 一 个 子 阵 ， 


0 
干 信号 的 舒缓 能 力 P? 


号 的 相干 性 影响 ， 
解 的 过 程 中 相 比 MD 算法 也 能 大 大 降低 
。 但是， 这 也 带 来 了 阵列 孔径 的 
同时 ， 当 存在 非 


声 时 ， 
， 将 会 导致 


党 声 子 空间 ， 从 而 得 出 信号 DOA， 这 就 是 
正 的 协 方差 矩阵 分 解 (smoothed matrix 


DOA 
接收 信号 


的 自 | 


写 空 间 


度 与 分 间 平 


向 平滑 信号 避 


其 中 : X*(1) 


™ 


代表 X00 的 复 
前 向 平滑 信和 号 


Y(1)=JX"(7) 


共和 圈 ， 
X() 的 协 方差 入 


J 和 矩阵 满足 J* 
E 车 即 为 及 =R,， 而 根据 式 


(7) 可 得 后 向 平滑 的 信号 协 方差 矩阵 为 
R, =E[Y(DY™(D)]=JRiJ=JA'R A J+JQY 


其 中 : 


ls 


则 第 m 个 阵 元 的 输出 信号 为 
x (1) = > ee| - jd, sin& j+g, (7?) (2) 
其 中 : si (7) 表示 接收 到 的 第 个 源 信号 ;9 为 第 个 信号 源 
的 入 币 ? 度 ; 8g,,(?) 示 表示 m 信道 的 非 均 匀 高 斯 噪声 。X 为 
号 波长 ， 定 义 9. 的 方向 向 量 为 
X=AS+G (3) 
其 中 : 
X=[(D, (人 We DY; 
A=[a(0),a(@),…,a(0.)]; 
S=[s (Ds (0) esa (OF ; (0 
G=[g,(),8,(1),…, gu (OD ; 

非 均 匀 品 声 集合 G~ CN(0,@ ; Q 为 非 均匀 噪声 的 协 方差 
功率 矩阵 ， 当 快 拍 数 为 了 时 根据 式 (3) 可 得 接收 信号 协 方差 矩 
阵 为 : 

R=E[X()X"()]=ARA'+Q (5) 


同时 ， 


R,=diag[P,P,,..., 


R] 


9 


让 


中 凡 = ELS()3"(O]= 了 SS" 表示 发 射 信号 的 协 方差 矩阵 ; 


RPR 为 第 1! 个 信号 功率 ， 


可 利 


根 据 


rank(Ri) =rank(R: -Qi) ， 其 
则 以 上 问题 可 以 ] 


则 新 构成 的 数据 协 方差 矩阵 
及 与 及 都 可 以 表示 
Ft 及 与 R, 分 别 
而 Qi 与 Q@ 代表 前 
接着 可 以 通过 凸 优化 原理 


尺 表 前 向 


F J 以 R, 
号 源 数 ( M >>K )， 
低 秩 矩 阵 ， 且 噪声 所 对 应 
| Ri 的 低 秩 
行 重 构 ， 从 而 抑制 非 
及 =A+Q: 


中 Q, 为 3 


特性 与 凸 优 
匀 匀 的 高 斯 噪 


Rup 变 为 


Rv =[R, 及 ] 


为 例 ， 


s 为 及 =R+Qi 与 局 =R,+Q, 的 


滑 算法 对 于 相 
1, OC-SMD 算法 是 将 整个 天 线 阵 列 都 看 
畜 。 则 前 向 待 平 滑 的 接收 信号 即 为 X(4) ， 而 后 
以 表示 为 


(9) 


=1。 


(10) 


(11) 


NN 


与 后 向 空间 的 平滑 无 噪声 
向 与 后 向 空间 的 非 均匀 
分 别 来 求 出 前 向 
假设 信号 
则 rank(R)=K <<M 
的 元 素 值 即 为 A 
化 原理 


与 


前 述 为 


Ey 
R,=R,+Q, 


th 


正定 矩阵 。 


噪声 协 方差 矩阵 。 
与 后 向 的 无 噪 
号 模型 中 的 阵 元 数 远 
， 所 以 Ri 是 
的 主 对 角 线 元 素 ， 
来 对 其 主 对 角 线 元 素 进 
体 算法 如 下 : 


协 方差 ， 


可 得 
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min rank{R. — Q,} 
st. R,-Q:>~0, (12) 
Qi eZ 


Z* 代表 正定 矩阵 集合 。 但 是 上 述 的 目标 函数 是 非 凸 的 ， 
所 以 没 办 法 在 如 此 庞大 的 NP-hard 问题 中 找到 全 局 最 优 解 C2 
31。 可 以 通过 核 范 数 来 凸 近似 矩阵 的 秩 ， 所 以 上 述 问题 可 以 
转换 为 


min IA. -Ql, 
st. R-Q.>~0, (13) 
Qe eZ 
而 矩阵 X 的 核 范 数 1XL 定义 为 
[Xl =tr(V XX) (14) 


所 以 它 又 可 以 等 价 为 矩阵 的 迹 ， 又 因为 
所 以 式 (13) 又 可 以 等 


Qu = diag(qs) (gst =[qit, ge», gu 了) 


价 为 : 
上 tr(R — Qi)= min [tr(Ri) -tr(Qi)]= 
max tr(Q)=max (Lg.,) (15) 
st. PR.—diag(q4)>0, 
qu >0 , (qs =[qirs gar. qr] ) 
式 (15) 中 L= 呈 4…,4lw ， 则 式 (13) 简 化 为 了 半 正 定 规划 问 


问题 (15) 可 以 利用 cvx29 等 工具 解 出 9* ， 继 而 可 求 得 前 
向 平滑 空间 的 无 噪声 信号 协 方差 矩阵 R 。 
与 前 向 平滑 空间 的 算法 一 样 。 通 过 凸 优化 处 理 可 以 算出 


分 别 得 到 Ri 与 R, 之 后 ,利用 式 (11) 即 可 得 到 去 掉 非 均匀 
高 斯 噪声 后 的 数据 协 方差 矩阵 Rw。， 然 后 对 该 修正 的 矩阵 进 
行 SVD， 就 可 得 到 修正 后 的 信号 子 空间 与 噪声 子 空间 ， 最 后 
通过 最 小 化 搜索 得 出 信号 源 的 方向 。 

根据 给 出 的 阵列 模型 与 MD 算法 的 子 阵 划分 方法 ， 则 数 
据 协 方差 矩阵 有 min(P,q) 个 相互 独立 的 矢量 ， 并 且 还 要 满足 
q>K 才能 实现 对 天 个 相干 信号 源 的 DOA 估计 。 即 MD 算法 
在 以 上 条 件 下 ， 最 多 可 估计 min(P,q) 个 相干 信号 源 。 而 
OC-SMD 算法 将 整个 线 列 看 为 一 个 子 阵 ， 协 方差 矩阵 具有 更 
多 的 独立 矢量 ， 最 多 可 估计 M -1 个 相干 信号 源 ， 解 决 了 MD 
算法 带 来 的 孔径 损失 的 问题 。 

算法 步骤 如 下 : 

a) 将 整个 天 线 阵列 看 做 是 一 个 子 阵列 ， 分 别 算出 前 向 平 
滑 空间 的 信号 协 方差 矩阵 与 后 向 平滑 空间 的 信号 协 方差 矩阵 
RR 与,; 

b) 利 用 凸 优 化 原理 分 别 来 求 出 前 向 与 后 向 的 无 噪声 信号 
协 方差 矩阵 Ri 与 R,; 

c) 利 用 SMD 算法 得 到 去 掉 非 均匀 高 斯 噪声 后 的 数据 协 
方差 矩阵 ,然后 对 该 修正 的 矩阵 进行 SVD 就 可 得 到 修正 后 的 
信号 子 空间 与 噪声 子 空间 。 

dj) 通过 最 小 化 搜索 实现 DOA 估计 。 
3 ”仿真 分 析 

为 了 验证 OC-SMD 算法 的 性 能 , 将 采用 等 间距 直线 阵列 
进行 仿真 实验 ， 入 射 信号 均 为 远 场 窗 带 信号 ， 信 源 数目 己 知 
目 均 为 相干 信号 。 将 会 与 MD 算法 ,Zi-SVD 算 法 以 及 SS-CRM 
算法 进行 各 个 方面 的 比较 。 
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3.1 相同 条 件 下 MD 和 OC-SMD 算法 检测 相干 信号 源 数目 

对 比 

在 阵 元 数 为 8 的 等 间距 直线 阵列 中 ， 信 品 比 为 15 dB， 
快 拍 数 为 1024。 当 子 阵 数 为 5， 每 个 子 阵 的 阵列 数 为 4 时， 
MD 算法 可 以 实现 最 多 相干 信号 源 的 估计 ， 此 时 可 以 估计 4 
个 相干 信号 源 。MD 算法 的 入 射 角 为 -40°、-20°20°.40? ， 
OC-SMD 算法 的 入 射 角 为 -60?、-40°、-20°.0°.20°?.40°.60? 。 图 2 是 
两 种 算法 在 八 阵 元 均匀 直线 阵 条 件 下 的 空间 谱 。 
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(b) OC-SMD 算法 在 八 阵 元 均匀 直线 阵列 中 的 空间 谱 
(b) Spatial Spectrum of OC-SMD algorithm in Eight-Array element 


uniform Linear Array 
2 MD 与 OC-SMD 算法 的 空间 谱 对 比 
Fig.2 Comparison of Spatial Spectra between MD and OC-SMD 


algorithms 
图 2 可 以 看 出 ， 利 用 OC-SMD 算法 ，8 元 阵 可 以 同时 
检测 7 个 相干 信号 源 ， 而 MD 算法 最 多 只 能 检测 4 个 相干 信 
号 源 。OC-SMD 算法 在 检测 相干 信号 源 数目 上 明显 多 于 MD 
算法 ， 很 好 地 克服 了 孔径 损失 的 问题 。 
3.2 不 同 信 噪 比 下 四 种 算法 的 空间 谱 对 比 

在 阵 元 数 为 8 的 等 间距 直线 阵列 中 , 有 两 个 相干 信号 源 ， 
入 射 角 分 别 为 -5° 和 0" 。 令 信 品 比分 别 为 15 dB、5 B 和 -5 dB， 
快 拍 数 为 1024，MD 算法 的 子 阵 数 为 5， 每 个 子 阵 的 阵列 数 
为 4。 图 2 是 四 种 算法 的 空间 谱 对 比 。 
由 图 3 可 得 以 下 几 点 分 析 结 果 : 
a) 由 图 3(a) 可 以 看 出 ， 当 信 噪 比 为 15 dB 时 ， 四 种 算法 
都 能 准确 地 分 辨 出 两 个 信 源 的 入 射 角 ; 
图 3(b) 可 以 看 出 ，5dB 这 样 较 小 的 信 噪 比 对 MD 以 
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及 LI-SVD 算法 的 影响 都 很 大 ， 从 其 空间 谱 图 像 中 已 经 不 能 


c) 


所 阁 ， 等 : 基于 修正 后 矩阵 分 解 的 最 优 协 方差 DOA 估计 


3.3 不 同 信 品 比 下 的 四 种 算法 的 均 方 根 误差 
在 阵 元 数 为 8 的 等 间距 直线 阵列 中 ,有 两 个 相干 信号 源 ， 


够 分 辨 出 角度 差 为 5° 的 两 个 信号 源 DOA。 较 小 的 信 噪 比 同 


算法 从 图 像 来 看 ， 依 然 可 以 分 辨 出 
出 了 更 好 的 分 辨 性 能 ; 


时 也 对 OC-SMD 和 SS-CRM 算法 产生 了 影响 ,但 是 , OC-SMD 


个 信号 源 的 方向 ， 表 现 


SS-CRM 算法 已 经 无 法 从 其 空间 谱 


的 入 射 角 
但 是 在 -5 dB 的 环境 下 ， 依 然 能 够 在 其 空间 谱 图 像 中 
个 相干 信 源 的 入 射 角度 。 所 以 , 在 低 信 噪 比 情况 下 ， 


影响 ， 


辨别 出 


度 , 而 OC-SMD 算法 虽然 


图 3(@) 可 以 看 出 ， 当 信 噪 t 


进一步 恶化 到 -5 dB 时 ， 


也 会 受到 低 信 噪 比 带 来 的 


OC-SMD 算法 的 分 辨 性 能 更 好 。 
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图 像 中 辨别 出 两 个 信号 源 
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(b) 信 噪 比 为 5 dB 时 的 空间 谱 


(b) Spatial Spectrum with a signal-to-noise ratio of 5 dB 
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Fig.3 Spatial spectral contrast 
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入 射 角 分 别 为 -5 和 0 。 令 快 拍 数 为 1024，MD 算法 的 子 阵 
数 为 5， 每 个 子 阵 的 阵列 数 为 4。 针 对 [-8:2:12] 范 围 内 不 同 的 
信 噪 比 , 通过 1000 次 的 蒙特 卡 洛 实验 , 得 出 四 种 算法 的 均 方 
根 误差 〈root mean square error, RMSE) 变 化 情况 。 


率 更 高 ， 但 是 低 信 噪 比 对 其 性 能 影响 还 是 很 大 ; 
c) SS-CRM 算法 虽然 在 低 信 噪 比 条 件 下 的 DOA 估计 误 
差 很 小 ， 但 是 低 信 噪 比 对 其 性 能 的 影响 也 不 可 忽略 ; 


估计 均 方 根 误差 /(*) 


0 2 4 
信 噪 比 /dB 


图 4 不 同 信 噪 比 条 件 下 的 均 方 根 误差 变化 情况 


Fig.4 Change of root mean Square error under different 


signal-noise-ratio conditions 


图 4 可 得 以 下 几 点 分 析 结 果 : 


日 


a) MD 算法 作为 传统 算法 ， 对 信 噪 比 的 改变 非常 敏感 ， 
在 低 信 品 比 环境 中 对 相干 信 源 的 入 射 角 估 计 的 误差 非常 大 ; 


b) L1-SVD 算法 虽然 较 MD 算法 ， 在 低 信 品 比 时 的 准确 


d) OC-SMD 算法 的 RMSE 在 不 同 信 噪 比 条 件 下 始终 是 


四 


种 算法 中 最 低 的 , 尤其 在 0 dB 以 下 的 低 信 噪 比 区 间 中 , 准 


确 


率 较 之 其 他 三 种 算法 更 高 。 


3.4 不 同 快 拍 数 下 四 种 算法 的 均 方 根 误差 

在 阵 元 数 为 8 的 等 间距 直线 阵列 中 ,有 两 个 相干 信号 源 ， 
入 射 角 分 别 为 -5° 和 0 。 令 信 噪 比 为 5 dB，MD 算法 的 子 阵 
数 为 5， 每 个 子 阵 的 阵列 数 为 4。 针 对 [100:100:1100] 范 围 内 


不 同 的 快 拍 数 , 通过 1000 次 的 蒙特 卡 洛 实验 , 得 出 


的 RMSE 变化 情况 。 
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图 5 不 同 快 拍 数 下 的 


四 种 算法 的 均 方 根 误 差 


Fig.5 Change of root mean square error under different snapshots 


日 


a) 随 着 快 拍 数 的 增加 ， 四 


图 5 可 得 以 下 几 点 分 析 结 只 


-人 小: 


个 算法 的 均 方 根 误差 都 在 逐渐 


减 小 。 尤 其 是 MD 算法 ， 其 准确 率 在 低 信 噪 比 条 件 下 并 不 理 


想 ， 而 它 对 于 快 拍 数 的 依赖 更 是 大 大 影响 了 角度 分 辨 性 能 ; 
b) Zi-SVD 算法 虽然 在 快 拍 数 降低 的 时 候 其 RMSE 较 
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MD 算法 更 低 , 但 是 因为 它 在 5 dB 这 种 低 信 噪 比 环境 下 的 误 ” 功 抑 制 非 均匀 噪声 对 DOA 算法 带 来 的 巨大 影响 。 实 验 结果 
差 本 来 就 较 大 ， 所 以 在 图 4 中，Li-SVD 的 DOA 估计 性 能 依 表明 ， 相 比 起 MD 算法 、Zi-SVD 算法 以 及 SS-CRM 算法 ， 
然 较 差 ; OC-SMD 在 低 信 噪 比 、 小 快 拍 数 条 件 下 的 DOA 估计 误差 更 
c) SS-CRM 算法 虽然 在 快 拍 数 变化 的 情况 下 依然 能 够 保 小 。 同 时， 在 非 均 匀 噪 声 环境 下 ， 拥 有 比 SS-CRM 算法 更 强 
持 较 为 稳定 的 RMSE， 但 是 从 图 中 可 以 直观 地 看 出 当 快 拍 数 ”的 噪声 抑制 能 力 。 
较 低 的 时 候 ，SS-CRM 算法 的 DOA 误差 还 是 比较 大 ; 参考 文献 : 
d) OC-SMD 算法 与 其 他 三 个 算法 相 比 , 其 RSME 在 小 快 一 
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